Conservation Properties in the Time-Dependent Hartree Fock Theory 
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We discuss the conservation of angular momentum in nuclear time-dependent Hartree-Fock cal- 
culations for a numerical representation of wave functions and potentials on a three-dimensional 
cartesian grid. Free rotation of a deformed nucleus performs extremely well even for relatively 
coarse spatial grids. Heavy ion collisions produce a highly excited compound system associated 
with substantial nucleon emission. These emitted nucleons reach the bounds of the numerical box 
which leads to a decrease of angular momentum. We discuss strategies to distinguish the physically 
justified loss from numerical artifacts. 

PACS numbers: 21.60.Jz, 25.70.-z 



Timc-dcpcndcnt Hartree- Fock (TDHF), originally pro- 
posed by Dirac has found widespread applications in 
various areas of physics due to the overwhelming devel- 
opment of computational power. It is employed, e.g., 
as the variant time-dependent density functional the- 
ory H in atomic, molecular and cluster physics, see 
e.g. [H, HJ. It has enjoyed application in nuclear dy- 
namics since more than thirty years [5| as a microscopic 
approach to various dynamical scenarios in the regime 
of large amplitude collective motion, like fusion exci- 
tation functions, fission, deep- inelastic scattering, and 
collective excitations; for early reviews see, e.g., p, jTj . 
With the steady upgrade of computational power, three- 
dimensional TDHF calculations without any symmetry 
restriction became possible and renewed the interest in 
nuclear TDHF as seen from an impressive series of re- 
cent publications @, |, M, EH, El El EI EE El E3 • A 
crucial aspect in nuclear TDHF is that nuclei are freely 
moving objects such that all conservation laws (energy, 
momentum, angular momentum) apply. Conservation of 
energy and momentum is a basic feature which has been 
tested for all existing codes. Conservation of angular mo- 
mentum has not yet been studied and that is the topic 
which we want to address in this paper. 

The calculations employ grids in coordinate space. 
Their finite spacing and box size destroy translational 
and rotational symmetry which, in turn, can spoil con- 
servation of momentum and angular momentum. The 
major destructive mechanism comes from matter which 
tries to leave the computational box but is hindered by 
the boundary condition. In the course of time develop- 
ment, higher energy components appear in the nucleon 
wave functions, representing outgoing "particles" from 
the nucleus. The further time development may there- 
fore be affected by their reflection from the boundary or, 
in the case of periodic boundary conditions, re-entry from 
the neighboring cells. That can change the total angular 
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FIG. 1: Simple illustation of the effect of boundary condi- 
tions on total angular momentum. The boxes indicate the 
computational boundaries and the central dot the reference 
point for the angular momentum. For periodic boundary con- 
ditions (left) it is easy to even revert the sign of the particle's 
angular momentum. For approximately reflecting boundary 
conditions (right) the situation is not quite as pronounced: 
for the case vi there is no change by reflection, while for «2 
the sign also changes. 



momentum, as is illustrated in Fig. [TJ Clearly the bound- 
ary can even change the sign of a particle's contribution 
to the angular momentum around the center of the cell. 
In this work we consider periodic boundary conditions. 
The case of reflecting boundaries behaves qualitatively 
similar. 

The static and time-dependent Hartree-Fock equations 
are solved on a cartesian three-dimensional mesh with- 
out any symmetry requirements. The grid spacing was 
1 fm, and the Skyrme energy functional [181 ] was em- 
ployed with the parametrization SLy6 [h1 ] (for the pur- 
poses of this work, the particular choice of Skyrme force 
is irrelevant). The minimum set of time-odd terms to 
assure Galilei invariance [H, HI Hl| was included. The 
spatial derivatives are calculated using the Fast Fourier 
Transform and periodic boundary conditions are em- 
ployed, except for the Coulomb potential, which is calcu- 
lated with boundary conditions at infinity as described in 
[22| . The time stepping employs a sixth-order Taylor ex- 
pansion of the time evolution operator U(t,t + At) = 
exp [— \h{t + At/2) At/h] , with the mean fields at the 
half step estimated by a third-order expansion using the 
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mean field at time t. 

This method of time development is non-unitary, so 
that the orthonormality of the single-particle wave func- 
tions is not guaranteed. Nevertheless, we find that the 
calculation is quite stable and accurate for several thou- 
sand time steps At sw 0.2 fm/c, in the sense that the 
particle number changes by less than 0.1%. The total 
energy also is conserved to a fraction of an MeV. When 
instability then sets in, there is a rapid drift of particle 
number and energy, so that the conservation properties 
are quite good checks for the accuracy and stability of 
the calculation. 

Test case will be a collision of two 16 O nuclei which in- 
volves truly large amplitudes and carries a large amount 
of excitation energy. As argued above, the boundaries of 
the numerical grid influence the conservation laws. In or- 
der to disentangle its effects, we consider three different 
setups: 

1. A small grid with 24 x 32 x 32 fm 3 . 

2. A doubled grid of 48 x 64 x 64 fm 3 , so that the colli- 
sion is surrounded more generously by empty space. 
Obscrvables are summed over only the smaller grid 
of case (H}, which allows to distinguish the exact 
physical loss from the artifacts of the smaller grid. 

3. Finally, the small grid plus an absorbing bound- 
ary conditions which are arranged within an ab- 
sorbing layer iV a b s = 6 cells wide in each direc- 
tion. In this layer, a mask function M(n x , n y , n z ) = 
M x (n x )M y (n y )M z (n z ) is applied to the wave 
functions after each time step where Mj(nj) = 

cos ((N=bs + l-ni)7r/2N abs )°' 25 , n; = l...N abs . 
See [23[ for details. 

Before analyzing rotational motion, we have checked, 
of course the conservation of energy and total momen- 
tum. Both quantities are conserved very well with rela- 
tive fluctuations staying at the order of 10 -4 . A detailed 
analysis of translational motion and of the physical in- 
terpretation of the TDHF single-particle energies can be 
found in Let us just mention here that for a nucleus 
moving freely in any direction on the grid the momentum 
is conserved to an accuracy of better than 10 -4 . 

For the case of angular momentum, the situation 
turned out to depend strongly on the excitation of the 
system. We therefore discuss two types of calculations: 
single cranked nucleus and heavy-ion collisions. 

A. Single cranked nucleus 

We produce a rotating nucleus by solving the cranked 
static Hartree-Fock equations 

(h-uj x ) <t>k(r) = ek<f>k{r), k = l...A (1) 

with u> the prescribed angular frequency of rotation 
about the x-axis (for simplicity we omit spin dependence, 
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FIG. 2: Relative deviations of the angular momentum expec- 
tation value (full curve) and the three principal moments of 
inertia (smallest: dotted, intermediate: dashed, and l largest: 
dot-dashed curve) from the temporal average during rotation 
of a 24 Mg nucleus. The abscissa denotes the rotation angle 
calculated from the instantaneous tensor of inertia. 



though spin is included in the calculations). The rotating 
states 
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then are exact solutions of the TDHF equations, so that 
the numerical solution should show simply a rotating nu- 
cleus with no extraneous motion added. 

It is clear that this is a much more demanding test for 
the numerical solution, since the cartesian grid is incom- 
patible with rotational motion and effectively the grid 
spacing expands and shrinks by a factor of y/2 as the 
nucleus rotates through 45 and then 90 degrees. 

As an example, we show here the deformed nucleus 
24 Mg cranked with uo = 2MeV/?i. In this case, rotation 
turns out to be almost completely of rigid-body type; the 
observed angular momentum of 7.54 h corresponds to a 
moment of inertia of about 98% of the rigid-body value. 
To judge the accuracy of the rotation in the numerical 
solution, we examine the expectation value of the angular 
momentum J x , which should be strictly conserved, and 
the principal moments of inertia, which in the absence of 
any internal excitation should also be constant. Fig. [2] 
shows the fluctuations of these quantities as functions of 
simulation time. 

The most striking result is the excellent quantitative 
description of rotation in spite of the coarse Cartesian 
grid, the variations being of the order of 10~ 4 . The 
angular-momentum expectation value clearly is corre- 
lated with the angle and shows regular variations with 
a 45 structure. The variations in the moments of inertia 
show a less regular pattern and appear to be influenced 
by the periods of the internal vibrations of the nucleus 
from which seemingly a tiny amount is exited during the 
rough sliding over the grid. 
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FIG. 3: Angular momentum as a function of time for collisions of 16 0+ ie O at E cm = 25MeV (left) and 125 MeV (right). The 
different cases are: "small": calculation in a grid of 24 x 32 x 32 fm 3 , "large": calculation in a grid doubled in total size in 
every direction, "restricted": same as large, but angular momentum is summed up only over the small grid; "absorb": small 
grid with a 6 fm absorbing layer around the boundary. 



B. Heavy-ion reactions 

The situation is quite different in the more violent case 
of heavy-ion collisions. The substantial excitation leads 
to emission of nucleons which, in turn, causes problems. 
As was shown in the introduction, particles crossing the 
boundaries can generate large spurious changes in the 
total angular momentum and it becomes quite difficult to 
separate the correct physical loss of angular momentum 
carried away by the emitted particles from the spurious 
numerical effect. The results presented in Fig.[3]illustrate 
the problems which are of quite general nature. Note 
that the angular momentum J x is perpendicular to the 
reaction plane. 

The initial condition consists of two ground-state 16 O 
nuclei with a cm. energy of either 25 or 125 MeV. All the 
initial angular momentum thus comes from the relative 
motion. For the higher energy the impact parameter was 
b = 4.8 fm, corresponding to J x as 33ft, while for the 
lower energy b = 2.8 fm corresponding to J x as 8.6/L 
These values were chosen to have the same distance of 
closest approach in the pure Coulomb trajectory. 

For both energies, the two nuclei stay fused at the end 
of the calculation. The boundary problems are more seri- 
ous for the higher energy because of the higher excitation 
leading to stronger emission of particles. It is apparent 
that in the larger grid the reduction in angular momen- 
tum starts about 100 fm/c later as compared to the small 
grid. The curve labeled "restricted" is computed in the 
large grid while the angular momentum is collected in the 
small grid. This should indicate the true loss of angular 
momentum from the small grid for the time span before 
emitted particles come back from the larger boundary. 
Clearly the "small" calculation has the largest loss and 
it becomes even unreasonable at about 250 fm/c. It is re- 
assuring that the curve for the absorbing boundary stays 
quite close to the "restricted" calculation and shows a 
reasonable monotonic decrease throughout. 

The total reduction of about 30% at a collision en- 
ergy of 125 MeV is surprisingly large in view of the fact 



that only 1.7 nucleons are absorbed. These nucleons thus 
carry a comparatively large share of angular momentum. 
At the lower energy of 25 MeV, the total change in angu- 
lar momentum is not as dramatic but by no means negli- 
gible, still exceeding 10% while 0.4 nucleons are emitted. 
Besides the quantitative difference, the general pattern 
are very similar. Again, the loss sets in later for the 
larger grid and the calculation with the absorbing layer 
appears to be a reasonable approximation to the "true" 
loss. 

In this paper, we have analyzed the conservation 
of total angular momentum in nuclear TDHF calcula- 
tions. The calculations used a coordinate-space represen- 
tation of wave functions and potential fields on a three- 
dimensional cartesian grid without any symmetry restric- 
tion. The full Skyrme interaction was taken into account. 
Conservation of energy and momentum was tested (but 
not detailed here) and found to be well matched within 
a relative error of only 10~ 4 . 

The results for angular momentum depend on the dy- 
namical scenario. Free rotation of a deformed nucleus is 
surprisingly well described. Although the cartesian grid 
spoils rotational symmetry, we find that the deformed nu- 
cleus rotates steadily over the grid of 1 fm spacing with 
variations of angular momentum and moments of inertia 
of the order of 10~ 4 . This is the same quality as found 
already for translational momentum and energy. 

The case of nucleus-nucleus collisions is less well- 
behaved. The compound system is heavily excited. This 
leads to substantial emission of nucleons which in the se- 
quel reach the bounds of the numerical box where reflec- 
tion or periodic copy (depending in the grid model) lead 
to a substantial reduction of angular momentum. Com- 
paring calculations on different grids (small box, large 
box, absorbing bounds), we have worked out that the 
loss is to a large extent physical because the emitted nu- 
cleons are very energetic and carry away a comparatively 
large amount of angular momentum. Artifacts from the 
boundary come into play as soon as nucleons travel back 
into the reaction zone. This happens the later the larger 
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the grid and it can be effectively avoided when using ab- 
sorbing boundary conditions. Both "solutions" , larger or 
absorbing grid are somewhat expensive. One may live 
with a small grid if one confines the analysis to the early 
time evolution of the compound system (to evaluate the 
doorway effects). In any case, the conserved quantities 
should be checked carefully for each new dynamical sce- 
nario. 

The results discussed here are, of course, only indica- 
tive and may vary quantitatively for other TDHF codes. 



The specific discretization of the equations of motion will 
affect the accuracy of describing an isolated rotating nu- 
cleus, while the treatment of the boundary conditions 
will strongly influence the boundary problems addressed 
above. 
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